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ABSTRACT 


A formula  is  derived  for  interpolation  between  output  samples  of  an  FFT 
i*e»,  in  the  frequency  domain.  Such  a formula  is  useful  for  obtaining 
greater  frequency  resolution  when  two  coarse  FFT  outputs  are  available. 
Consideration  is  also  given  to  the  effect  of  such  interpolation  on  a 
weighted  FFT. 
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In  digital  signal  processing  applications  involving  spectral  measure- 
ments, frequency  resolution  is  limited  by  the  length  of  the  input  sample  stream 
and  by  available  data  processing  capacity.  However,  in  applications  such  as 
doppler  radar,  it  often  happens  that  high  frequency  resolution  is  needed  only 
over  regions  which  are  small  compared  to  the  sampling  rate,  whereas  coarser 
information  will  do  for  the  rest  of  the  spectrum.  Such  increased  resolution 
may  be  achieved  in  a number  of  different  ways: 


One  may  increase  the  order  of  the  filter  bank  Fast  Fourier  Transform  (FFT) 
Normally  the  order  is  chosen  to  be  a radix  2 number  in  order  to  optimize 
the  efficiency  of  the  FFT.  Thus,  if  the  order  is  doubled,  the  processing 
time  more  than  doubles  before  any  information  becomes  available. 


If  the  frequency  regions  of  interest  are  preselectable,  high  resolution 
may  be  obtained  if  they  are  small  (see,  for  example,  p.  390  of  £lj  ).  The 
breakeven  point  turns  out  to  be  a region  size  of  half  the  sampling  rate. 
This  method  is  not  universally  applicable  in  that  no  coarse  information 
is  made  available  for  the  full  range  of  frequencies,  the  frequency  regions 
swat  be  selected  apriori  (rather  than  as  a result  of  processing  coarse 
data),  and  the  region  sizes  must  be  simple  fractions  of  the  sampling  rate 
such  as  1/3,  1/12,  etc. 


(c)  Another  alternative  is  to  form  two  low  order  (say  N)  FFTs  from  the  input 
data,  and  then  combining  (over  the  frequency  regions  of  interest  only)  to 
obtain  the  higher  resolution  2Nth  order  FFT.  The  first  FFT  would  use  the 
even  numbered  samples,  and  the  second  the  odd,  as  in  the  mechanization  of 
a decimation- in- time  FFT.  The  drawback  with  this  method  is  that  the  out- 
puts of  the  primary  FFTs  will  furnish  coarse  spectral  data  aliased  2:1. 
Thus,  with  the  first  half  of  the  spectrum  folded  into  the  second  half, 
the  primary  FFTs  will  yield  intelligible  coarse  information  only  if  the 
spectral  regions  of  interest  are  folded  into  regions  which  are  clear  (i.e 
which  contain  noise  only).  Even  when  this  happens  to  be  the  case,  the 
aliasing  reduces  the  coarse  data  signal-to-noise  ratio  by  3 db.  Again, 
no  coarse  spectral  data  would  be  available  before  2n,  data  points  were 
collected  and  processed.  | — r |_|  * ^ | ft  * g 


Finally,  one  could  perform  an  Nth  order  FFT  on  the  first  N samples  and 
another  on  the  next  N.  The  outputs  of  either  of  these  can,  of  course, 
be  used  as  coarse  data  to  determine  regions  where  greater  resolution  is 
required.  (In  doppler  radar  applications,  for  example,  targets  could  be 
detected  with  the  coarse  data).  The  question  now  arises;  can  the  two  Nth 
order  FFT  outputs  be  easily  combined  to  produce  greater  resolution  in 
frequency  areas  of  interest?  The  answer  is  "yes".  We  now  derive  the 
formula  to  be  used  in  this  data  combination. 


Let  the  two  sequential  FFT  output  vectors  be  FFTl  and  FFT2  respectively 
Let  the  2n  inputs  be 


Let  the  high'  resolution  spectrum  be 


Then,  thinking  in  terms  of  a decimation- in- frequency  FFT,  we  see  that  the 
even  outputs  are 


where 


are  the  even  members  of  the  output 
The  odd  members  are 
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and  thus 
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The  frequency  domain  interpolation  formula  we  require  consists  of  eqs.  (6)  and 
(12).  The  formula  is  exact,  but  in  order  to  gain  the  computational  advantage 
which  is  desired  in  interpolation  formulae,  one  would  want  to  sum  over  only  a 
limited  range,  say,  for  where  L,  and  L^are  close  to  l<  . Note 

that  the  importance  of  the  £ to  the  spectral  estimate  as  approximately 

inversely  proportional  to  Jt-k  . In  fact,  if  N is  large,  we  see  that 
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whenever  we  use  only  a few  neighboring  d^  such  that 


If  this  is  done,  we  may  set 
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Eq.  (14)  may  be  used  as  an  approximation  to  the  odd  numbered  values,  eq.  (12), 
However,  for  a given  application  there  is  not  a particularly  great  advantage 
in  doing  so.  The  coefficients  of  <5.^  given  in  eq.  (12)  are  easily  obtained  on 
a computer. 

At  this  point  one  might  ask  how  big  the  interval  should 

be,  and  if  the  terms  outside  this  interval  may  indeed  reasonably  be  neglected. 
Since  the  sum  in  eq.(12),  or  equivalently,  in  the  approximating  eq.(14),  is 
finite,  there  is  no  question  of  divergence.  However,  the  basic  question,  "How 
good  is  the  approximation?",  is  probably  best  answered  by  deriving  a formula 
for  the  finite  impulse  response  (FIR)  filter  weights  corresponding  to  the  inter- 
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polated  values,  and  Chen  checking  their  frequence  response. 
Using  nearest  neighbors  only,  eq.  (12)  becomes 
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where 
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where 
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Viewing  eq.  (17)  as  a convolution  sum,  we  see  that  the  FIR  weights  are 

|(3*  = s^w 

Generalizing  the  definition  (16)  to 
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we  can  write  an  expression  for  the  FIR  weights  using  the  nearest  M neighbors. 
We  first  set 
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Comparing  eqs.  (20)  with  eqs.  (15)  and  following,  we  easily  deduce  that  the 
weights  are 
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If  the  approximation  leading  to  eq.  (14)  is  used,  we  obtain 
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for  the  FIR  weights  using  the  M nearest  neighbors.  It  is  desirable  to 
normalize  the  weights  as  derived  in  eqs.  (19),  (22),  or  (23)  in  order  to 
achieve  a specific  gain  at  the  center  of  the  filter.  In  the  following 
example  we  show  magnitude  responses  of  a few  candidate  interpolating  filters. 
They  are  all  normalized  to  OdB. 


Example 


Suppose  N = 16,  2n  * 32,  and  k = 2.  I.e. , we  desire  the  spectral  estimate 

based  on  two  sequential  16-point  FFT  output  vectors. 

(Such  a small  order  is  chosen  here  for  illustrative  purposes  only.)  Figure  1 
shows  how  the  high  resolution  spectrum  is  formed  using  nearest  neighbors  only. 
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Figure  1.  Interpolation  with  nearest  Neighbors  only. 

Figure  2 shows  the  result  of  using  nearest  neighbors  only,  formula  (19). 

In  this  case  N Dolph-Chebyshev  weights  with  40  dB  sidelobe  suppression  were 
used  [2].  Cali  these  weights  { UJW  } K=  O,  - - M - 1 } • 

The  equivalent  set  of  2N  weights .becomes 
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The  filters  corresponding  to  the  neighboring  even-numbered  spectral  estimates, 
G-tf  and  , are  represented'  by  dashed  curves  in  the  figure.  As  expected, 
using  N weights  leads  to  rather  poor  filters  after  combining  FFTl  and  FFT2. 
Note  that  such  preweighting,  which  is  commonly  used  for  better  filter  shaping, 
is  easily  incorporated  into  the  formulas  derived  above.  Since  the  weights  go 
with  the  input  samples,  equations  (6)  and  (12)  remain  unchanged,  and  (22) 
becomes 
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The  magnitude  response  of  these  W is  what  is  graphed  in  figure  2.  This 
kind  of  weighting  is  appropriate  for  the  two  N-point  FFTs,  but  is  undesirable  for  a 
2N-point  FFT.  The  appropriate  weighting  would  consist  of  2n  weights.  Figure  3 
shows  Clif  , and  in  this  case.  The  interpolating  filter's  sidelobe 

structure  is  now  much  more  appealing,  and,  of  course,  the  even  numbered  filters 
are  of  just  the  desired  Dolph-Chebyshev  type.  Figure  4 shows  the  same  situation, 
but  in  this  case  the  four  nearest  neighbors  on  each  side  have  been  used  (formula 
(24)  with  jurjj  as  2n  Dolph-Chebyshev  weights.)  As  expected  the  result  is 
much  improved.  The  drawback  with  using  2n  weights  is  obvious:  each  of  the 
primary  FFTs  receives  a lopsided  set  of  weights.  However,  since  these  coarse 
spectra  are  used  for  quick-look  purposes  only,  their  degraded  quality  may  often 
be  tolerable.  Figure  5 shows  the  magnitude  response  of  FFTlQ,  the  first  filter. 

The  dashed  curve  shows  the  filter  shape  with  N Dolph-Chebyshev  weights. 

Conclusion 

This  paper  has  derived  a set  of  equations,  (6)  and  (12  or  14),  useful  for 
interpolating  between  weighted  (or  unweighted)  output  samples  of  sequential  FFTs, 
i.e.,  in  the  frequency  domain.  The  quality  of  the  interpolation  has  been  illus- 
trated by  means  of  a numerical  example,  which  also  showed  the  consequences  of  two 
different  weighting  approaches.  It  is  challenging  to  consider  that  a compromise 
may  exist  between  these  two  approaches,  l.e.»  a weighting  sequence  which  would 
Improve  both  the  high  resolution  and  the  quick- look  spectra  as  compared  to  the 
results  of  either  of  the  above  approaches. 


Figure  :3.  (Magnitude  Response  for  the  CX^; Filter  using 
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Figure  5.  Magnitude  Response  for^N  lopsided  Weights 
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